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HIGHLIGHTS 


► Torrefaction based co-firing system was proposed. 

^ The kinetics of torrefied biomass has been investigated experimentally. 

^ The performances of co-firing with varying biomass substitutions have been studied. 
^ No obvious reduction of boiler load was found when 100% coal switched to biomass. 
^ Net C0 2 and NO* emissions reduced with increasing of biomass substitutions. 
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Torrefied biomass has several benefits, such as higher energy density, good grindability, higher flowabil- 
ity and uniformity. The process of torrefaction moves the chemical and physical properties of raw bio¬ 
mass close to that of bituminous coal, which allows co-utilization with high substitution ratios of 
biomass in the existing coal-fired boilers without major modifications. In this study, a torrefaction based 
co-firing system was proposed and studied. Devolatilization and char oxidize kinetics of the torrefied bio¬ 
mass have been investigated experimentally. CFD modeling of co-firing with varying substitutions of 
torrefied biomass in a pulverized coal boiler have been carried out. To figure out the boiler performance 
when co-firing torrefied biomass, five different cases were involved and simulated, coal only, 25% bio¬ 
mass, 50% biomass, 75% biomass, and 100% biomass on thermal basis, respectively. The results showed 
torrefaction is able to provide a technical option for high substitution ratios of biomass in the co-firing 
system. The case-study pulverized coal boiler could be fired 100% torrefied biomass without obvious 
decreasing of the boiler efficiency and fluctuation of boiler load. More positively, the net C0 2 and the 
NO* emissions significantly reduced with increasing of biomass substitutions in the co-firing system. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Co-firing biomass with coal in traditional coal-fired boilers for 
electricity production represents a cost effective and efficient 
renewable option that promises reduction in effective C0 2 emis¬ 
sions, reduction in NO* emissions, and several societal benefits 
[1,2]. Many different types of biomass can be co-fired with coal, 
such as wood, residues from forestry and related industries, agri¬ 
cultural residues, and various biomasses in refined form such as 
pellets. Biomass is able to penetrate all energy sector markets, 
but technical and economic constraints still limit its general 
deployment. Currently, technical challenges associated with 
co-firing are: (1) cost-effective methods for getting more and wider 
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varieties of fuels into the boiler, (2) occurrence of insufficient gas 
mixing or stratified flows in the boiler, (3) fouling and corrosion 
of the boiler, (4) continuation of fly ash utilization, and (5) impacts 
on performance of the flue gas cleaning. The economics of biomass 
co-firing are significantly impacted by the twin factors of deregula¬ 
tion and environmental requirement. Those non-technical barriers 
are as follows: (1) lacking of financial incentives, and uncertain fuel 
price/availability, (2) legislative aspects (utilization of fly ash in ce¬ 
ment, determining green share of electricity produced, and unclean 
emission legislation), (3) public perception of co-firing of biomass/ 
waste, and (4) getting the permits through. Uncertainties in both 
arenas make investment in most, or all, biomass co-firing technol¬ 
ogies difficult to justify. Therefore, it is expected an efficient but 
sustainable biomass thermal conversion schemes. Further on, 
new EU regulations force power sectors to use more biomass, such 
as Large Combustion Plant Directive (LCPD), National Emission 
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Nomenclature 



A 

pre-exponential factor for char oxidation (s -1 ) 

H 

hydrogen content in solid fuel (wt.%, daf) 

At 

pre-exponential factor for devolatilization (s _1 ) 

m p 

particle mass (kg) 

c 

carbon content in solid fuel (wt.%, daf) 

m P , o 

initial particle mass (kg) 

C D 

drag force coefficient 

O 

oxygen content in solid fuel (wt.%, daf) 

Csw 

swelling factor 

Re 

Reynolds number 

d p 

current particle diameter (m) 

T m 

mean particle and gas temperature (K) 

dp,o 

initial particle diameter (m) 

Vi 

volatile yield factor 

Do 

diffusion rate coefficient (m 2 /s) 



E c 

the activation energy of char oxidation (kj/mol) 

Greek symbols 

Ei 

activation energy of devolatilization(kj/mol) 

Pp 

particle density (kg/m 3 ) 

fv, 0 

the fraction of volatiles in the particle 

P 

fluid phase density (kg/m 3 ) 

fw, 0 

the fraction of moisture in the particle 

Pp 

particle velocity (m/s) 

g 

gravity acceleration constant, 9.8 m/s 2 

P 

fluid phase velocity (m/s) 


Ceilings Directive (NECD), and Integrated Pollution Prevention and 
Control Directive (IPPCD) [3]. Those regulations further promote 
large percentage biomass will be used by power sectors. 

The co-firing technologies employed in power sector are 
broadly classified into three types: direct co-firing, indirect co-fir¬ 
ing, and gasification based co-firing [4,5]. In direct co-firing system, 
the biomass is feed into the mills to be firstly pulverized as bio¬ 
mass powder in the mill system, and then introduced into boiler 
either separately or blend with coal [6-8]. Indirect co-firing in¬ 
volves the installation of a completely separate biomass boiler to 
produce low-grade steam for utilization in the coal-fired power 
plant prior to being upgraded, resulting in higher conversion effi¬ 
ciencies [4]. Gasification based co-firing system involves the gasifi¬ 
cation of solid biomass and combustion of the product fuel gas in 
the furnace of the coal-fired boiler, and offers a high degree of fuel 
flexibility [9-11 ]. Currently, direct co-firing is the most popular op¬ 
tion for biomass applications in industry, due to the relatively low 
investment cost of turning an existing coal power plant to a co-fir¬ 
ing plant [4,12]. However, the majority of direct co-firing installa¬ 
tions were operated with less than 10% (thermal basis) of biomass 
co-firing ratios, since the technical and economic constraints [4,5]. 

2. Torrefaction based co-firing system 

2A. Concept of torrefaction based co-firing system 

In order to increase the possibility of co-firing application of 
biomass, the torrefaction as a pre-treatment technology for co-fir¬ 
ing was reported by the Energy Research Centre of the Netherlands 
(ECN) in 2005, which makes biomass is able to compete with fossil 
fuels, and thus provides an option for the direct co-firing with a 
large percentage of biomass without major threats [13]. 

Torrefaction is a thermo-chemical process conducted at 200- 
300 °C, typically, 1-h residence time, during which biomass par¬ 
tially decomposes giving off volatiles and producing the remaining 
solid as a final product [12,14]. Loss of the tenacious nature of raw 
biomass is mainly coupled to the breakdown of hemicellulose ma¬ 
trix, which bonds the cellulose fibers in the raw biomass, and de¬ 
crease the length of these fibers during the depolymerization 
process [15]. Generally, torrefied biomass has several benefits 
[13-18]. First, torrefied biomass contains, typically, 70% of its ini¬ 
tial weight and 90% of the original energy content [19]. As the re¬ 
sult, the energy density of biomass increases about 1.3 times after 
torrefaction process, so that the stability of the boiler loads can be 
achieved, and the transportation cost of torrefied biomass will be 
also reduced about 1.3 times. Moreover, the higher fibrous in 
raw biomass makes the grinding of the raw sample difficult and 


energy consuming [14]. Compared with the raw biomass, torrefied 
biomass can be characterized by the anhydrous weight loss, and 
the grinding energy greatly decreases with weight losing [20]. 
For example, a reduction of 85% of power consumption was maxi¬ 
mally observed for torrefied willow [21 ]. Thirdly, in contrast to the 
raw biomass, torrefied wood particles have higher flowability and 
good fluidization behavior, which is an important performance for 
an easy injection of biomass powder into boiler, because the torr¬ 
efied biomass has a brittle behavior with decreased strength, smal¬ 
ler particles sizes, and lower moisture content. Finally, torrefied 
biomass has an increased uniformity in product quality. For in¬ 
stance, it has been found that woodcuttings, demolition wood, 
waste wood have quite similar physical and chemical properties 
after torrefaction, so that it is possible to use a wide range of fuels 
and reduces season influences [13]. Accordingly, the promising 
chemical and physical properties of torrefied biomass provide a 
co-utilization option with high ratios of biomass in existing coal- 
fired boilers without major modifications. However, the combus¬ 
tion behavior, operation performance and emissions of co-firing 
system based on biomass torrefaction are not clear. 

In this study, the torrefaction based co-firing system in pulver¬ 
ized coal boilers has been proposed with the aim of 100% fuel 
switching, as schemed in Fig. 1. Raw biomass was pretreated by tor- 
refaction process previously and the torrefied biomass was grinded 
in the existing but separate coal-mills to produce the biomass pow¬ 
ders, which were finally introduced into combustion chamber by 
primary air. Torrefied biomass powders were then combusted sep¬ 
arately in the upper burners when in co-firing cases. The bottom 
burners are designed either for coal particles or biomass particles, 
which can be switched to use biomass only in the pure biomass 
combustion cases, as shown in Fig. 1. To figure out the performance 
of co-firing torrefied biomass in a front-wall pulverized coal boiler, 
five different cases were involved and simulated, coal only, 25% bio¬ 
mass, 50% biomass, 75% biomass and 100% biomass cases respec¬ 
tively. In which, the substitution ratios of biomass are on thermal 
basis, and the 100% coal case has been selected as the base case 
for validation. More detailed investigations covering the profiles 
of temperature, major gas components, pollutants and unburned 
carbons have been done. The results of numerical analysis with 
aim of 100% fuel switching were discussed in this paper. 

2.2. Fuel properties of torrefied biomass 

In this work, torrefaction process has been done previously to 
get torrefied biomass, which derived from its raw biomass with 
20% of volatile content releases on dry basis. Proximate and ulti¬ 
mate analysis data of raw biomass, torrefied biomass and coal 
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Fig. 1. Scheme of torrefaction based co-firing system. 


Table 1 

Analysis data of raw biomass, torrefied biomass and coal particles. 


1.6 ■ 



Units 

Raw biomass 

Torrefied biomass 

Coal 


1.4 - 

Moisture 

4.2 

4.2 

5.6 

1.43 



Proximate 

analysis (dry basis) 





1.2- 

Ash 

1.95 wt.% 

1.95 

7.31 

25.06 

G' 


Volatile 

79.72 wt.% 

79.72 

61.87 

27.85 

rz 

1.0 - 

Char 

18.33 wt.% 

18.33 

30.82 

47.09 

'w' 


LHV 

18.0 MJ/kg 

18.0 

20.5 

24.0 

U 

se 

- r 

00 

o 

Ultimate analysis (dry basis) 




(j 


C 

50.32 wt.% 

50.32 

55.43 

60.26 

1 


H 

6.02 wt.% 

6.02 

4.90 

3.97 


0.6 - 

O 

41.42 wt.% 

41.42 

31.83 

8.33 

<3 


N 

0.26 wt.% 

0.26 

0.50 

1.27 


0.4 - 

S 

0.03 wt.% 

0.03 

0.04 

1.11 




particles are listed in Table 1. The proximate analysis confirms that 
the torrefaction dramatically reduced the volatile to fixed carbon 
ratio and improves the energy density. Moisture content of the 
raw material is lower than the torrefied material, which could be 
caused by the different origin of the material. 

The Van Krevelen diagram gives the relation between O/C and 
H/C for different solid fuels, as shown in Fig. 2. In which, the torr¬ 
efied biomass and coal were studied in this work; the bituminous, 


0.0 - 


Raw biomass 

♦ 


Coal 

(this work) 


Sub-bituminous 


0.2 - 


0.0 


1 

•• 

.t 

Bituminous Torrefied biomass 
(this work) 


0.2 0.4 0.6 

Atomic O/C (daf) 


0.8 


Fig. 2. Van Krevelen diagram for different solid fuel (Note: (1) Raw biomass is the 
parent species of torrefied biomass; (2) The reference data of bituminous and sub- 
bituminous are from IFRF data bank). 


sub-bituminous and raw biomasses are used for reference. From 
Fig. 2, obviously, the raw biomass has lower energy density than 
coal. However, promisingly, it shows that torrefaction changes 
the elementary composition of biomass and moves it towards coal. 
Besides, breakdown of hemicellulose matrix with weakening of the 
tenacious nature of the biomass, which bonds the cellulose fibers 
in biomass, and decrease the length of these fibers during depoly¬ 
merization process. Therefore, torrefied biomass is able to be grin¬ 
ded by the existing coal milling systems without modifications. 

Co-firing a biofuel in existing pulverized-coal fired boilers sets 
high demands for grinding of biomass [6]. The critical factor is 
the particle size of the biomass, which is necessary to ensure prop¬ 
er combustion efficiency and a stable flame. In this work, after tor- 
refaction, it is assumed that biomass was grinded separately by an 
existing coal-mill. Then the biomass particles and coal particles 
sieved and selected were in same range of 62-150 pm. 


2.3. Experimental studies of torrefied biomass 
2.3. t. TG and DTG tests 

To further identify the advances of chemical properties of torr¬ 
efied biomass, thermogravimetric tests were performed under in¬ 
ert or atmosphere conditions. In this work, the operating 
procedure includes, dry process, which is a heating process with 
constant heating rate (20°C/min) from room temperature to 
105 °C under nitrogen, and then maintains the isothermal period 
at 105 °C about 10 min; pyrolysis process, which is also a heating 
process with constant heating rate (20°C/min) from 105 to 
900 °C under nitrogen, and then maintains the isothermal period 
at 900 °C about 10 min; char oxidized process, which starts nitro¬ 
gen cooling period from 900 to 800 °C under nitrogen gas, and then 
switches from nitrogen to air, and maintains the isothermal period 
for 10 min under air. 
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2.3.2. Kinetics tests 

To obtain more realistic kinetics data, devolatilization tests at 
higher heating rate were carried out in an isothermal plug flow 
reactor (IPFR) [22]. The test conditions were set similar to industrial 
pulverized fuel applications, such as high temperature and high 
heating rate. The reactor is possible used to represent conditions 
found in full-scale combustion processes with heating rates on 
the order of 10 4 -10 5 °C/s, maximum temperature of 1400 °C, and 
particle residence time between 5 and 1500 ms. A flow of 100 g/h 
of biomass powder was guaranteed to maintain stationary condi¬ 
tions for a time test of at least 30 min. The fuel conversion was eval¬ 
uated by comparing the composition of the collected solid residue 
with that of the parent fuel by means of the ash tracer method. 

After devolatilization under the nominal temperature of 900 °C 
and residence time of 300 ms conditions that assure the complete 
devolatilization but prevent the graphitization of the carbonaceous 
matrix, the char of torrefied biomass obtained have been tested in 
the ThermoGravimetric Balance to evaluate the char oxidation 
kinetics. Small char samples (3 mg) have been tested in order to 
limit material and thermal resistances during the tests. Three-time 
tests with a constant heating rate (20 °C/min) were carried out 
[ 22 ]. 

3. Description of boiler and operating conditions, and modeling 

3.1. Boiler introduction 

The case-study boiler, OP-650, is a front-wall pulverized coal 
boiler with maximum electricity production capacity of 
220 MWe. In this work, the full boiler geometry was modeled, 
and the cross sectional area of domain was 15.598 x 8.958 m 2 
and the height of the domain was 39.453 m, as shown in Fig. 3. 
The mesh system consisted of 796,372 structured hexahedral cells. 
The mesh was refined at the location of the burners, where the 
ignition and most of the combustion reactions occurred. The 
super-heaters on the top of boiler were also considered. There 
are 24 burners were arranged on front wall, all the burners as 4 
rows and 6 burners equally spaced in each row. The amount of pri¬ 
mary air controlled the rate of combustion in the chamber. The 
amount of secondary air controlled combustion efficiency. Suffi¬ 
cient OFA was added for the complete oxidation of any unburned 
combustibles and OFA ports were located on the rear wall, there 
were 10 injection nozzles. In addition, 8 vertical slots-air (front 
wall) and 2 rows protection-air (rear wall) were used for protecting 
the membrane walls against over temperature operation and influ¬ 
ence of combustion products. 

3.2. Operating conditions 

In this work, five different cases were simulated and studied. 
One case of pure coal, three cases of co-firing with different substi¬ 
tution portions of biomass (25%, 50%, and 75% on thermal basis), 
and one case of pure biomass were modeled and studied. The oper¬ 
ation data of 200 MWe pure coal case was collected in the power 


Table 2 

The operation parameters of boiler at 200 MWe. 

Electrical power of block (MWe) 

Boiler efficiency (%) 

Pressure in the boiler drum (MPa) 
Temperature in the boiler drum (°C) 
Temperature of the feed water (°C) 

Pressure of the feed water (MPa) 

Capacity of the feed water (t/h) 


plant, as shown in Table 2. To obtain the operating parameters of 
three co-firing cases and the pure biomass case, several assump¬ 
tions have been made based on the pure coal case: (1) total energy 
capacities are kept as 200 MWe, (2) total energy inputs of fuel are 
513 MW/h, (3) the boiler operated with excess air ratios of 1.15, (4) 
all 24 burners are operated with average energy capacity and same 
conditions in all co-firing cases and the pure biomass case, (5) the 
portions of PA, SA, OFA, and protection air are kept same as the 
pure coal case. The operating parameters of all five cases are listed 
in Table 3. 

3.3. Basic mathematical models 

To obtain as precise as possible numerical analysis results, Sol- 
idworks software was used for the air supply system simulation to 
predict the velocity vectors and the temperature of combustion air. 
The details of air supply system simulations have been reported by 
Li et al. [23]. All the velocity vectors of combustion air obtained 
from the air system simulation were inputted as boundary condi¬ 
tions in Fluent by user-defined methods. 

After the air supply system simulations, Ansys Fluent 13.0, was 
used to simulate combustion, fluid and particle flow, as well as 
heat and mass transfer inside the furnace [24]. The standard k-s 
model with standard wall function was used for calculating the ef¬ 
fect of turbulent flow inside furnace [25]. Lagrangian particle tra¬ 
jectories of the pulverized coal particles were calculated 
throughout the computational domain. The discrete ordinates 
(DOs) radiation model was used to simulate radiation heat trans¬ 
fer. Simulating coal combustion mainly involved the modeling of 
devolatilization, the homogeneous combustion of gaseous vola¬ 
tiles, and the heterogeneous oxidation of solid char particles [26]. 
The dry process, which is also important, starts with a progressive 
heating of the coal particle by convection and radiation. Wet com¬ 
bustion model was used for modeling the moisture release process 
during the heating stage [24]. 

Devolatilization is modeled to predict the volatile matter yield 
by the two-competing-rates model proposed by Kobayashi et al. 
[27] and the gas released during this process contributes to the gas¬ 
eous reaction. The char combustion is computed according to the 
kinetics/diffusion-limited model [28,29], where the surface reaction 
rate is determined either by kinetics or by a diffusion rate. In prac¬ 
tice, the operating temperature of furnace is high (above 1000 °C), 
so that volatile combustion happens at high Damkohler numbers 
(Da » 1). Accordingly, the global reaction limited by diffusion was 
used for combustion of gaseous volatiles, and calculated by Eddy- 
Dissipation model [24]. For the torrefied biomass particles, the data 
for the kinetic rate of devolatilization and char oxidation were mea¬ 
sured and presented in Table 4. For the coal particles, the kinetic 
parameters of devolatilization were obtained from Li et al. [30]; 
the kinetics of coal char oxidation were calculated based on Hurt 
and Mitchell’s correlations [31]. Besides, diffusion coefficients for 
all cases were calculated based on D 0 = 2.53 x 10“ 7 T^ 75 /d p [32]. 

3.4. NO x formation modeling 

Thermal NO* is formed from oxidation of atmospheric nitrogen 
at relatively high temperatures in fuel-lean environments and has 
strong temperature dependence. The formation of thermal NO is 
determined by extended Zeldovich mechanism [33]. Fuel NO* is 
formed from nitrogen bound in solid fuel, and the part of fuel- 
nitrogen is usually released as HCN and NH 3 accompanying with 
volatile releasing, which finally will be either oxidized into NO or 
reduced into N 2 [34]. For the coal, it is assumed that 90% of the 
nitrogen from the volatile will be converted to HCN, and the rest 
will form NH 3 [35]. While for biomass, oppositely, 90% of HCN 
and 10% of NH 3 were derived from the volatile [36]. Besides, in this 
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Fig. 3. Global geometry and grid structure of burners, OFA ports and vertical slots (dimensions in mm). 


Table 3 

The operating parameters of all five cases. 



Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

Combustion type 

Pure coal 

Co-firing 

Co-firing 

Co-firing 

Pure biomass 

Biomass substitutions (%, thermal basis) 

0 

25 

50 

75 

100 

Biomass burners 

/ 

B41-B46 

B31-B36 

B21-B26 

All 




B41-B46 

B31-B36 






B41-B46 


Total burners for biomass 

0 

6 

12 

18 

24 

Coal feed rates (t/h) 

79.2 

59.4 

39.6 

19.8 

0 

Biomass feed rates (t/h) 

0 

22.5 

45.0 

67.5 

90.1 

Primary air (N m 3 /h/°C) 

201170/100 

199726/100 

198017/100 

196309/100 

194600/100 

Secondary air (N m 3 /h/°C) 

224810/289 

223196/289 

221287/289 

219377/289 

217468/289 

OFA (N m 3 /h/°C) 

67865/289 

67377/289 

66801/289 

66225/289 

65648/289 

Protection air (N m 3 /h/°C) 

26301/289 

26112/289 

25889/289 

25665/289 

25442/289 

Vertical slot air (N m 3 /h/°C) 

33105/289 

32868/289 

32586/289 

32305/289 

30024/289 


study, it is assumed that the nitrogen contained in the char will be 
directly converted to NO, since its complex conversion mecha¬ 
nisms. Finally, prompt NO also was considered, which is formed 
by the reaction of atmospheric nitrogen with hydrocarbon radicals 
in fuel-rich regions of flames [37]. 

3.5. Discrete phases modeling (DPM) 

Both biomass particles and coal particles are modeled by the 
DPM approach. A key consideration in modeling trajectories of so¬ 
lid fuel is the various forces acting upon the particles. These initial 
simulations included specific terms for gravitational, buoyancy and 
drag forces. The trajectories, heat and mass transfer from/to these 
discrete particles are coupled with the continuous phase. The tra- 


Table 4 

Kinetic parameters for torrefied biomass and coal particles. 

Torrefied biomass 3 Coal b 

Devolatilization kinetics (at high heating rates) 


Al (s- 1 ) 

5.48E+03 

3.7E+5 

A2 (s” 1 ) 

2.15E+07 

1.46E+13 

El (kj/mol) 

59.8 

73.7 

E2 (kj/mol) 

147 

251 

VI 

0.62 

0.3 

V2 

0.70 

0.8 

Char oxidation kinetics 

A(s -1 ) 

56.3 

0.78 

Ea (kj/mol) 

78 

94.6 


a The kinetics of torrefied biomass were measured in this work. 
b The kinetics of coal were from [30,31]. 
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Fig. 4. TG, DTG and temperature profiles at heating rate of 20 °C/min. 


jectory of the discrete phase particle is calculated by integrating 
the force balance on the particle as motion equations [24]: 


du p _ 1 Sju C D Re 


dt 


pp d P 


(U - Up) + 


s(p p - p) 

Pv 


3.5A. Drag force coefficient, C D 

Based on the above equation, obviously, drag force is important 
for calculating the trajectory of discrete particle. And the drag force 
coefficient is a function of shape factor of particle [38]. The shape 
factor 0 , a typical indicator of irregularity of solid particles, was 
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Fig. 5. Profiles of temperature and gas components along the height of boiler. 
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Table 5 

Comparison of heat absorptions in different cases. 



Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

Biomass substitutions (%) 

0 

25 

50 

75 

100 

Heat absorbed by membrane walls (MW) 

229.73 

222.29 

215.73 

206.28 

191.47 

Heat absorbed by super-heaters (MW) 

38.19 

44.28 

47.15 

49.05 

50.95 

Total heat absorptions (MW) 

267.92 

266.57 

262.88 

255.33 

242.42 

Discrepancies (%) 

Base case 

-0.50 

-1.88 

-4.69 

-9.52 


§ 

'o* 


o 

u 



0 25% 50% 75% 100% 


Biomass substitutions (Thermal basis) 


4. Results and discussions 

4.1. Comparison of torrefied biomass and raw biomass 

For comparison, the TG, DTG tests of both raw biomass and torr¬ 
efied biomass at heating rate of 20 °C/min have been carried out, as 
shown in Fig. 4. Generally, the first peak of devolatilization is 
attributed to the hemicellulose decomposition, which is the more 
reactive chemical compound of biomasses. In Fig. 4, the earlier 
sub-peak in the DTG curves is at 215 °C for the raw biomass and 
around 280 °C for the torrefied biomass. The delay of this sub-peak 
observed in the torrefied materials could be explained in a signifi¬ 
cant modification of hemicelluloses structure during the torrefac- 
tion process. The main peak of devolatilization was observed 
around 360 °C both for raw biomass and torrefied biomass, which 
was attributed to the decomposition of cellulose. The decomposi¬ 
tion of lignin occurred in a wider range of temperature from 
250 °C to the end of the process. Besides, a higher DTG sub-peak 
of oxidation in torrefied biomass than that in raw biomass shows 
that torrefaction improved the reactivity of char. 


Fig. 6. C0 2 emissions at different substitution ratios of biomass. 


4.2. Kinetic parameters of torrefied biomass 


defined as a ratio of surface area of a sphere area with the same 
volume as the particle and the actual surface area. The shape of 
the biomass particles is affected by the fiber content, although 
the torrefied biomass has less fiber content, but its surface area 
is still higher than that of pulverized coal. The higher surface area 
enhances the drag force imposed onto biomass particles and the 
contact possibility of 0 2 and particle surface. For coal particles, it 
is widely agreed that the drag force was calculated by spherical 
drag law [24,39]. In this work, a non-spherical shape was assumed 
for torrefied biomass particles with a shape factor of 0.54, and the 
drag coefficient was calculated based on Haider and Levenspiel’s 
equation [38]. 


3.5.2. Swelling coefficient , C sw 

In terms of above equation, the drag force is affected by particle 
diameter (dp). The diameter of a particle changes during devolatil¬ 
ization according to swelling coefficient. The relationship of parti¬ 
cle diameter and swelling coefficient is expressed as [24]: 


d P (t) 

dp, o 


1 + (Ssw - 1) 


(1 /w, o trip) 

fv,o (1 — fw,o)rtip,o 


where ( \ is the ratio of the mass that has been devolatil- 

ized to the total volatile mass of the particle. The effect of swelling 
on drag force during the biomass devolatilization process should be 
taken into account, due to the significant fraction of volatile matter 
in biomass. In this work, swelling coefficient of torrefied biomass 
particles was also studied experimentally, which was calculated 
as the ratio of the average diameter of char particles to the average 
diameter of the initial biomass particles. Finally, the measured 
swelling coefficient of torrefied biomass is 1.8. According to the 
experimental studies of Fletcher [40], the swelling coefficient of 
pulverized coal was selected as 1.1. 


Two competing reactions are used to simulate the devolatiliza¬ 
tion step and predicted the high temperature volatile matter yield, 
as above mentioned. Briefly, the first reaction (Al, El, and VI) are 
used to prevail at lower temperature, while the second reaction 
(A2, E2, and V2) plays dominant role at higher temperature. Gener¬ 
ally, the volatile matter yield at high temperature is higher than 
the standard volatile matter revealed. The prediction of the volatile 
matter yield as function of the temperature is thus important for 
practice combustion models. The detailed kinetic data for biomass 
devolatilization are presented in Table 4. Once the pre-exponential 
factor and activation energy are known, the model provides the 
rate of volatile evolution and the related char and volatile yields. 
Finally, the average values of torrefied biomass char oxidation 
are also listed in Table 4. 

4.3. Validation 

In this paper, the predicted results of the pure coal case have 
been compared to the measured values. It is necessary to point 
out there are some differences between measured values and pre¬ 
dicted values due to different positions of them, the measured data 
were only available before air-preheater, but the predicted results 
were the Mass-Weight-Average values at the furnace exit. In terms 
of measured data, average concentrations of 0 2 , CO, and NO x were 
3.13%, 1124 mg/Nm 3 , and 535 mg/N m 3 (at 6% 0 2 ), respectively. 
The predicted results show good agreements with measured val¬ 
ues, and the predicted average concentrations of 0 2 , and CO were 
3.3%, 856 mg/N m 3 , and 530 mg/N m 3 (at 6% 0 2 ), respectively. In 
the real condition, the temperature of flue gas exiting the furnace 
is still high enough for oxidation reactions. Therefore, a lower 0 2 
content before the air-preheater than at the furnace exit would 
be reasonable. However, the predicted value of CO concentration 
was lower than the average of the measured values, probably 
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because the rate of oxidization of carbon in fly ash was faster than 
the CO consumption rate after the furnace exit. Therefore, it could 
be concluded that the adopted numerical models in present calcu¬ 
lations are reasonable for this study. 

4.4. Effects of biomass substitutions on temperature and gas 
components 

Profiles of temperature and main gas components along height 
of boiler with varying biomass substitutions are shown in Fig. 5. 
The boiler height used in Fig. 5 starts from the bottom of boiler do¬ 
main, as shown in Fig. 3. The bottom burners are located on height 
level of 8.3 m, the top burners are located on height level of 17.1 m, 
the OFA ports are installed on height level of 19.6 m, and the fur¬ 
nace exit is on height level of 27.6 m. 

The comparison of the flue gas temperature is clearly shown 
along the height of boiler with five different cases in Fig. 5a. 
Obviously, the flue gas temperatures distributions in those cases 


are almost shown in the same trends, which might be concluded 
that chemical properties of torrefied biomass and coal are close, 
and torrefied biomass can be a good option for replacing coal. From 
Fig. 5a, the average temperature of the pure coal case was obvi¬ 
ously higher than that of co-firing cases in the main combustion 
chamber (from the bottom burners to the top burners), while after 
top burners(especially, up to the OFA ports), the average tempera¬ 
tures of co-firing cases were higher than temperature of the coal 
case. The probably reason is, the motion process was faster than 
combustion process, due to the enhanced drag force, part of com¬ 
bustible components of biomass was conveyed upwardly to this 
zone by gas flow. Therefore, the unburned combustibles were oxi¬ 
dized in this zone, which produced a slightly higher temperature in 
co-firing cases. 

The 0 2 of air injected into furnace is quickly consumed with the 
combustion reactions happened, since the temperature inside 
furnace is enough higher to oxidize the volatile and char of solid 
fuels. In consequence, the 0 2 concentration decreased rapidly after 


Char contents (kg Char/kg Particle remains) 



(a) Injected from top burners 
(Case 5 /biomass only) 



i 

V • 


(c) Injected from 3 rd row burners 
(Case 5 /biomass only) 



(b) Injected from top burners 
(Case 1/coal only) 



(d) Injected from 3 rd row burners 
(Case 1/coal only) 


Fig. 7. Comparison of char fates and particles trajecoties of biomass case and coal cases. 
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Fig. 8. NO emissions in different substitution ratios of biomass. 

injection, as shown in Fig. 5b. In Fig. 5b, the 0 2 contents in all co¬ 
firing cases were obviously higher in the whole furnace, except the 
zone close OFA ports. The higher ratios of biomass, the oxygen dis¬ 
tribution was more similar with that of the pure coal case. The 
probable reason is the part of oxygen for biomass combustion 
was from oxygen content of biomass, so the oxygen devolatiliza¬ 
tion process from biomass partly slowed down the combustion 
process of biomass. Compared to coal char, the biomass char was 
more reactive, but obviously, the content of biomass char was low¬ 
er. Besides, turbulence flow pattern inside furnace might have 
some influences on oxygen concentrations in all co-firing cases. 

Actually, there are two routes of CO formation, either the oxida¬ 
tion of volatile or char oxidation. The CO concentration on the fur¬ 
nace exit is important to the combustion efficiency of furnace, and 
a lower CO concentration represents good combustion situation of 
boiler. The OFA have a significant influence on the CO concentra¬ 
tion, which was injected into furnace above burners for complete 
combustion. Thus, the CO concentration on the furnace exit can 
be used to judge if the amount of the OFA is enough for complete 
combustion. The CO concentrations in all cases along the boiler 
height are shown in Fig. 5c. In Fig. 5c, in all cases, the CO concen¬ 
trations increased with increasing of biomass substitutions, and a 
significant reduction was observed after the OFA ports. Interest¬ 
ingly, the CO concentrations in the main furnace zone (from the 
bottom burners to the top burners) of three co-firing cases and 
the pure biomass case were higher than that of the pure coal case. 
Therefore, it can be concluded that the major combustion zone of 
biomass was moved upwardly to the furnace exit. The reason 
might be the oxygen devolatilization process from biomass partly 
slowed down its complete combustion process, because of an 
important part of oxygen was from biomass itself. The other expla¬ 
nation might be the insufficient gas mixing or stratified flows post¬ 
poned a part of volatile combustion reactions. 

Finally, the characteristics of NO emissions were investigated 
for all cases. Total NO distributions along the boiler height in all 
cases are shown in Fig. 5d. It is found that the significant reduc¬ 
tions of total NO when switching the coal as torrefied biomass. 
Especially, the reduction of NO emissions is more clearly when 
the 50% of biomass on thermal basis was used in co-firing system. 
It seems that the thermal NO mostly formed in the upper zone of 
boiler since the slightly higher temperature when considered the 
temperature distribution in Fig. 5a, but simultaneously, it should 
be acknowledged that the short residence time of flue gas in the 


upper zone limited the formation of thermal NO, because the tem¬ 
perature of fuel gas was then quickly cooled down by the super¬ 
heaters. Moreover, the overall fuel-N contents decreased with 
increasing of biomass substitutions, due to the lower nitrogen con¬ 
tent of biomass compared to coal. Additionally, the fuel nitrogen in 
biomass volatile was mostly converted to NH 3 , which has a lower 
conversion ratio of NO when compared to HCN [7]. 

4.5. Effects of biomass substitutions on boiler efficiency 

Boiler load reduces when coal-fired boiler operated as co-firing 
mode, due to biofuel has a high water content [4]. When directly 
using a high moisture biofuel, the overall energy efficiency of the 
co-firing plant would be reduced unexpectedly. The process of tor- 
refaction moves the chemical and physical properties of the raw 
biomass close to that of bituminous coal with releasing almost 
all the moisture and less part of volatiles, which is expected to real¬ 
ize high biomass co-firing ratios without reducing energy effi¬ 
ciency. Therefore, the thermal performance of a torrefaction 
based co-firing boiler has been investigated in this work. 

Generally, a lower boiler load reduces the heat absorption of in¬ 
stalled heat exchangers, and therefore the temperature and pres¬ 
sure of main steam. In this work, to analysis the fluctuations of 
the boiler loads when the torrefied biomass replaced coal with 
varying biomass substitutions, the heat absorbed by the membrane 
walls and the super-heaters in different cases have been investi¬ 
gated and compared, as given in Table 5. The absorbed heat is de¬ 
fined as the integrations of heat flux over the whole area of studied 
walls/surfaces. In all cases, constant temperatures were kept as 
thermal boundary conditions of both the membrane walls and 
the super-heaters. From Table 5, it is observed that a relative less 
heat absorbed by membrane walls and a relative more heat ab¬ 
sorbed by super-heaters at a higher biomass substitution. This 
might be caused by the major combustion zone of biomass was 
moved upwardly to furnace exit, as discussed before. The discrep¬ 
ancies of total heat absorptions are also listed in Table 5. The larg¬ 
est discrepancy of -9.52% was found in the pure biomass case 
when compared to the pure coal case, which means that no more 
than 10% of boiler load was reduced if 100% coal was replaced by 
the torrefied biomass. Consequently, torrefied biomass provides a 
promising co-firing option with a higher boiler efficiency and min¬ 
or fluctuation of the boiler load. 

4.6. Effects of biomass substitutions on C0 2 emissions 

Biomass is an important renewable resource and the C0 2 re¬ 
leased during biomass combustion will be re-captured by the re¬ 
growth of the biomass through photosynthesis [41,42]. C0 2 emis¬ 
sion from the combustion of biomass fuel source is not assumed 
to increase the net atmospheric C0 2 levels [42]. Therefore, biomass 
is widely accepted as a potential solid fuel to reduce the net emis¬ 
sion of C0 2 . Typically, a large subcritical coal-fired utility plant pro¬ 
duce 900 kgC0 2 /MW h, and 750kgC0 2 /MWh for a supercritical 
plant [43]. To understand the different fates of C0 2 , capturable 
C0 2 and net C0 2 with unit of kg/MW h were compared based on 
the predicted data, as shown in Fig. 6. It is defined that the captur¬ 
able C0 2 derives from biomass combustion and the net C0 2 origins 
from coal combustion. Fig. 6, as expected, shows that the net C0 2 
reduces with increasing biomass percentages, while the capturable 
C0 2 increases simultaneously. There are 915 kgC0 2 /MWh recap¬ 
tured when 100% coal substituted by the torrefied biomass. When 
50% of biomass used in the co-firing boiler, the emission of net C0 2 
is 403 kg/MW h, which is less than half of emission of the net C0 2 
in the pure coal case. 

Moreover, the total C0 2 emission of the pure biomass case and 
three co-firing cases are slightly higher than that of the pure coal 
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case. Firstly, because the total energy inputs of fuel are same for all 
cases (see 3.2). According to the proximate and ultimate analyses 
of coal and torrefied biomass, it can be theoretically calculated that 
the total carbon content of fuel in the pure biomass case is slightly 
higher than that in the pure coal case. Secondly, the unburned car¬ 
bon contents in biomass ash were less than that of coal ash, since 
the biomass char is more reactive compared to coal char. It is com¬ 
monly acknowledged that the surface combustion of biomass char 
is much faster than of pulverized coal, due to enhanced mass trans¬ 
fer and intrinsic reaction kinetics. According to the prediction, the 
char burnout ratios of biomass in all co-firing cases were 100% and 
the total char burnout ratio for the pure biomass case was 99.95%, 
while the total char burnout ratio of coal for the pure coal case was 
97.75%. For clarify, the char contents of biomass and coal particles 
and their trajectories are clearly presented in Fig. 7. The biomass 
particles were almost combusted completely before leaving the 
domain when injected from top two rows burners; while coal 
particles showed similar trajectories, but the unburned carbon in 
coal ash is higher than that of biomass ash. 

NO molar reaction rates (mol/m 3 -s) 


4.7. Effects of biomass substitutions on NO emissions 

This substitution of biomass also has the potential to reduce the 
emission of air pollutants. Particularly, the emissions of NO x , which 
mainly due to the lower nitrogen content in the biomass. Fig. 8 shows 
the NO emissions with varying substitution ratios of biomass. Obvi¬ 
ously, the emissions of the fuel NO showed a linear and decreased 
with increasing of substitutions of biomass, and ratios of thermal 
NO/fuel NO were almost same in all cases, except the sudden 
decreasing of thermal NO at 50% biomass case. Therefore, in general, 
the total NO emission reduced with biomass substitution increasing. 
The total NO emission based on power output was 2.34 kg/MW h in 
the pure coal case, and almost half amount of which was reduced 
when 100% coal replaced by the torrefied biomass. 

To further understand the emission characteristics of the NO 
with varying biomass percentages in co-firing system, the forma¬ 
tion rate of thermal NO and fuel NO were studied separately, as 
shown in Fig. 9. It is observed that both the fuel NO and the ther¬ 
mal NO were mostly formed within an envelope of flames, since 
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Fig. 9. Comparison of NO molar reaction rates of biomass case and coal cases. 
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their formation rates primarily depend on the temperature and 
fuel/oxygen ratio. Near the burner region, due to the temperature 
of solid fuel rose up suddenly, the NH 3 or HCN were formed rapidly 
with releasing of the volatiles. Both in the pure biomass and the 
pure coal cases, the volumes with thermal NO formation rates 
above 0.0001 mol/m 3 -s are observed in the upper-center of the fur¬ 
nace when compared with that of fuel NO formation rates, as 
shown in Fig. 9, which means that an important part of the thermal 
NO was formed in the zone of above burners, where the interme¬ 
diates of the volatile-N subsequently reacted to form NO* under 
an oxidizing environment, because the OFA was injected proving 
the completely combustion air, but resulted in a lower emission 
of fuel NO. The size of the zone for the thermal NO formation rate 
was smaller compared to that of the fuel NO* formation for the 
pure coal case, which was similar with numerical results of Choi 
and Kim [44], this might because of the fuel NO formed more read¬ 
ily than the thermal NO during the coal combustion process. How¬ 
ever, compared to the pure coal case, thermal NO/fuel NO ratios 
increased with gradually increasing of the coal replaced by the 
torrefied biomass, which can also be observed from Fig. 8. Further¬ 
more, the major formation rates of the fuel NO in the pure biomass 
case showed a good corresponding to the envelope of flame, while 
the formation rate of the fuel NO in the pure coal case showed less 
order, as shown in Fig. 9c and d. 

5. Conclusions 

A torrefaction based co-firing system was proposed and studied 
systematically in this work. Torrefied biomass has several benefits, 
such as higher energy density, good grindability, higher flowability 
and uniformity in product quality. The process of torrefaction 
moves the chemical and physical properties of raw biomass close 
to that of bituminous coal, which allows co-utilization with high 
substitution ratios of biomass in the existing coal-fired boilers with¬ 
out major modifications. CFD modeling of torrefaction based co-fir¬ 
ing system in a pulverized coal boiler has been carried out with the 
aim of 100% fuel switching. To figure out the performance of co-fir¬ 
ing of torrefied biomass in a front-wall pulverized coal boiler, five 
different cases were involved and simulated, coal only, 25% biomass, 
50% biomass, 75% biomass, and 100% biomass, respectively. 

According to the results, torrefaction is bale to provide a techni¬ 
cal option for high substitutions of biomass in the co-firing system. 
The case-study boiler can be operated with a higher substitution ra¬ 
tio of substitution of the torrefied biomass without obvious 
decreasing of energy efficiency and fluctuation of boiler load. More 
positively, the net C0 2 and the NO* emissions significantly reduced 
with increasing of biomass substitutions. However, torrefaction 
based the co-firing system should be further studied from perspec¬ 
tive of other common issues, such as fouling, corrosion of the boiler, 
fly ash utilization, impacts on performance of the flue gas cleaning, 
and so on. In addition, the torrefaction based co-firing system also 
should be evaluated according to economic criteria for application. 
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